Elastic theory of quantum Hall smectics: effects of disorder 
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We study the effect of disorder on quantum Hall smectics within the framework of an elastic 
theory. Based on a renormalization group calculation, we derive detailed results for the degrees of 
translational and orientational order of the stripe pattern at zero temperature and carefully map 
out the disorder and length-scale regimes in which the system effectively exhibits smectic, nematic, 
or isotropic behavior. We show that disorder always leads to a finite density of free dislocations and 
estimate the scale on which they begin to appear. 

PACS numbers: 73.40.Hm, 73.20.Mf, 46.65.+g 



Recent experiments m discovered a pronounced 
anisotropy in the resistivity of the quantum Hall system 
near half filling of higher Landau levels (Landau level 
filling factor v > 3), suggesting that the ground state at 
these filling factors breaks orientational symmetry, ft is 
believed that this is a manifestation of the formation of 
striped quantum Hall states whose existence had previ- 
ously been predicted based on Hartree-Fock (HF) calcu- 
lations (for a recent review, see ||). The existence 
of striped states in higher Landau levels has also been 
supported by exact numerical diagonalization studies 
and by the experimental confirmation jlj of some explicit 
predictions [^],^| for their transport properties. 

The HF calculations predict that the ground state near 
half filling of higher Landau levels is a unidirectional 
charge density wave (CDW) whose period a is of the order 
of the cyclotron radius 0. In addition to orientational 
order, this state also breaks translational symmetry in 
one direction. For this reason, it is often referred to as 
a quantum Hall smectic state. Going beyond mean-field 
level by including quantum and thermal fluctuations, one 
expects a variety of phases with different degrees of trans- 
lational and orientational order 



Here we study the effect of disorder on quantum Hall 
smectics at zero temperature within the framework of an 
elastic theory. As a function of disorder strength, we 
map out in detail various regimes in which the system 
effectively exhibits smectic (with both translational and 
orientational order), nematic (with orientational order 
only), or isotropic behavior. A concise summary of our 
results is provided in Fig. 1. 

At low energies, the relevant gapless excitations are 
fluctuations of the displacement field u(r, t) associated 
with translational symmetry breaking and defined 
as the phase of the charge-density modulation p(r, t) = 
Po cos[fc(r|| + w(r,r))]. The position vector r = (rii,rj_) 
is represented by its components parallel and perpendic- 
ular to the CDW wave vector k with k = 2n/a. We note 
that as opposed to the stripe position described by it, the 
stripe width is a massive mode that can be neglected. 
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FIG. 1. Illustration of regimes in dependence of dis- 
order strength Ay and transverse length scale L± for 
q s ~ O.Olfc as is typical for experiments. In the smectic regime 



Lx « £5 



translational and orientational order are well 



pronounced, / <C k~ 2 and j< 1. In the weak-disorder case, 
Ay -C Ay, this is followed by a quasi-smectic regime with 
quasi-long-range translational order. In the strong-disorder 
case, Ay ^> Ay, a nematic regime follows, where transla- 
tional order is effectively short ranged. The strict validity of 
our analysis based on elastic theory ends at L±_ = L^ 100 , 
where free dislocations appear. Nevertheless, since in the ab- 
sence of dislocations orientational order would be long-ranged 
(i.e. £° rlont = oo; except for extremely strong disorder), ne- 
matic order possibly persists even for L± > Ll isloc . dc trans- 
port measurements on macroscopic samples probe the regime 
in the indicated area. 

The chiral nature of quantum Hall systems leads to a 
dynamic behavior of quantum Hall smectics which differs 
in important ways from that of ordinary smectic liquid 
crystals. The elastic action for quantum Hall smectics 
takes the form || (/3 denotes the inverse temperature) 

P u q ^ J ° 

where M is a dynamic coefficient and E(t) denotes the 
elastic free energy. The inverse powers of the transverse 
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wave vector q± in the dynamic term are specific to quan- 
tum Hall smectics. Most importantly, they imply that 
smectic order is stable against (weak) quantum fluctua- 
tions at zero temperature in the quantum Hall system, 
even in the absence of an alignment force orienting the 
stripes along a preferred direction [||. 

The elastic free energy is of the smectic form 0] 



(Z n ) = 



E(r) 



Bit Vnu - 



;(V^) 2 



-K{V\uf + £l(V ± u) 2 



(2) 



The constant C is determined by the requirement 
(Viiit) = which fixes the CDW wave length. The com- 
pression modulus B\\ and the bending modulus K are 
related by B\\ = 4fc 2 K @. For B± = 0, this free en- 
ergy is rotation invariant due to the anharmonic terms 
||. We have included a crystal field with tilt modulus 
B± which accounts for the presence of an alignment force 
orienting the stripe pattern along a preferred direction. 
Experiments suggest that this alignment force is quite 
weak, B± -C B» . The presence of the tilt modulus leads 
to an important (transverse) length scale in the problem, 
defined by q s = (B^/K) 1 / 2 . For q±^ ~> q s , the bending 
energy of the stripes dominates over the tilt energy while 
for q± <C q s , the tilt energy dominates. 

Various kinds of disorder may affect the quantum Hall 
stripes. A disorder potential U(r) is included via an 
impurity action S- lmp y = J U{v)p(v,r). It can be 

rewritten as S imp>v = (p /2) f rT [V(r)e iku ^ + c.c], 
where V is effectively locally correlated, (V(r)V* (r')) ~ 
WS(r — r'), even for long-range-correlated impurity po- 
tentials U(r) because of the oscillatory nature of the 
charge density p(r,r). The effective disorder strength 
W is given by W ~ (U(— k)[/(k)}. For quantitative esti- 
mates, we employ the impurity potential created by the 
ionized impurities (of areal density ni mp ) in the spacer 
layer which dominates the transport properties of two- 
dimensional electron systems (2DES). For uncorrelated 
impurity positions, the disorder potential experienced by 
the electrons in the 2DES has W — e 2 ni mp v 2 ,(k). Here, 
Vd(q) = 2ire~ qd / e(q)q is the screened Coulomb interac- 
tion, e(q) the dielectric function, and d the distance of 
the impurity layer from the 2DES. 

While the bare action contains only the point disor- 
der U(r), the action of the renormalization group (RG) 
generates a second type of disorder. This effective disor- 
der corresponds to a random tilting force acting on the 
stripe pattern and is described by the additional term 



J Y\i du a] ex P i - E ^iK*] 
a ^ a 

+ 4f[ TT X [ A "( V ^ r ' ^)X V ^(r, r')) 



+Ay cos[fc(u Q (r, t) - up(r, t'))] 



(3) 



Here, Ay = (\/2)p 2 ] W measures the strength of the point 
disorder while Ah denotes the strength of the tilt disor- 
der. 

Since translational order is stable against quantum 
fluctuations at zero temperature |5), one might attempt 
to treat the disorder potential perturbatively by expand- 
ing the action in Eq. (|^) to quadratic order in u, setting 
Ah = 0, and neglecting dislocations. In the resulting 
quadratic theory it is straight forward to calculate the 
correlation functions 



0(r) 



(Kr')-«(r' + r)] 2 ), (4a) 
([VWr')-V>(r' + r)] 2 ), (4b) 



S imp ,h = f r T /i(r)V_Lu(r, t) in the action with a local Qv = [Ay/A^] 1 ' 5 <? s , The disorder strength separating 



which describe the degrees of translational and orienta- 
tional order, respectively. Neglecting quantum fluctua- 
tions which are additive in perturbation theory, one finds 
at zero temperature /(ry) = \ J q± (k 2 A v /B 2 Qp[l — 
(1 + Q ||r|| ) e -Quni] and f(r x ) = \ £ (fc 2 Ay/i? 2 Qf )[1 - 
cos(gj_r^)]. Here r\\ (r±) are vectors in the paral- 
lel (perpendicular) direction and Q||(gj_) = [{B±_q\ + 
Kq^/B^ 1 ' 2 . The analo gous expressions for g(r) differ 
from those for /(r) merely by an additional factor q\ in 
the integrand. 

One readily finds that /(r) is infrared divergent for 
any nonzero r, irrespectively of whether the tilt modulus 
B± vanishes or not. Because of the multistability of the 
disorder potential, perturbation theory breaks down on 
scales where / > a 2 [fTof , i.e. on all finite scales. Below, 
we perform a RG calculation which properly accounts for 
the multistability by a significant downward renormaliza- 
tion of the disorder strength on large scales. 

The perturbative expression can nevertheless be used 
to estimate the characteristic transverse scale gy, below 
which the renormalization sets in. To do so, we de- 
mand f(qy 1 e±) ~ a 2 with the provision that the inte- 
gral over q± in the perturbative expression for / be self- 
consistently restricted to q± > qy . The explicit expres- 
sion for qy depends on the disorder strength. In the weak- 
disorder regime qy <C q B , one finds qy = [Ay/Ay] 1 / 2 g s , 
while in the strong-disorder regime qy 3> q s one obtains 

V6, 



correlator {h(r)h(r')) = Afj<5(r — r'). 

We proceed in a standard fashion by averaging the free 
energy over the disorder potential using the replica trick 
which leads to the impurity-averaged replicated partition 
function 



these two limits is given by A c v = AttB^B^/^K 

We now turn to the RG calculation, restricting our- 
selves to zero temperature where the smectic state is sta- 
ble in the pure system. As the approach neglects topo- 
logical defects such as dislocations, the results are valid 
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only up to the scale L^ sloc wnere dislocations start to 
appear. This length will be estimated below. Our anal- 
ysis builds on previous RG studies for classical smectics 
P,pl| , adapting them to two dimensions where disorder is 
known to play a peculiar role p"2[ . We employ a conven- 
tional coarse-graining procedure, where all displacement 
modes in the range A < \q±\ < k with A = ke~ l are 
integrated out. (A length rescaling is omitted.) Treat- 
ing both the point disorder and the anharmonic terms 
in the elastic action to one-loop order leads to the flow 
equations 



dAy 
dl 

dA h 

dl 
dB K 



= -aiAy, 
= a 2 A v + a 3 A h , 
= -3a 3 £||, 



dl 
dK 

~~dT = 2a ^ 



(5a) 
(5b) 
(5c) 
(5d) 



where a\ (A) 



/c 4 AA v /47rB2Q3(A) and a 2 (A) 



c(fc 2 /A 2 )ai(A) with c a constant of order unity. More- 



B}/ 2 A h /16irK 5 / 2 A 3 for A > q s and a 3 







over, a 3 (A) - ^ 
for A <^ q B . The couplings ai and a 2 arise from the point 
disorder while a 3 is due to the anharmonic terms in the 
elastic action. We note that the tilt modulus B± remains 
unrenormalized. The renormalization of -By and K is due 
to the combined effect of disorder and anharmonic elas- 
ticity. The characteristic disorder scale qy can be recov- 
ered from the flow equations by the condition ai(gy) = 1 
which defines the effective onset of the renormalization 
of Ay. We first solve the flow equations ignoring the 
anharmonic terms in the elastic action (i.e. a 3 = 0). As 
discussed below, these terms should lead only to minor 
modifications. Solving the flow equation we find that 
the effective disorder strength Ay ~ A 2 flows to zero for 
A <C qv, while the tilt disorder strength A^ increases 
logarithmically. Due to this decay of Ay the coupling 
ai remains finite for vanishing A, which indicates that 
the renormalization of disorder is captured consistently 
on one-loop level. 

The detailed results depend on the disorder strength. 
In the weak- disorder regime qv <C q s or Ay -C Ay, we 
find for the translational order in the soft (perpendicular) 
direction 



(q v r ± ) 2 \n(q s /qy), 



-1 



/( r ^) ~ p { (lvr±_) 2 ln(l/gyr ± ), q' 1 < r ± < q v \ 



ln 2 (gyrj_ 



r± > qy 1 



(6) 



As opposed to the perturbative result, / is now finite 
for any finite r± and we can identify a correlation length 
^trans f Qr translational order by /(C^ ans ) = a 2 and 
find ^ t j'' ans ~ q^ 1 . The slow In 2 increase of / in the regime 



r± £Y ans , typical of two-dimensional disordered sys- 
tems [ }L2| , implies that the system exhibits "quasi" - 
long-range translational order, with power-law Bragg 
peaks with logarithmically varying powers. The orien- 
tational correlation function g remains finite, g(r±) ~ 
(qy/k) 2 \n(q s /qy) <C 1, for arbitrarily large lengths r± » 
qy 1 . Hence the correlation length £™ lent for orientational 
order, defined by <7(£° rient ) = 1, is infinitely large. 

For sufficiently strong disorder, we enter the strong- 
disorder regime qy 3> q s or Ay ^> Ay where 




(qv/q s ) 3 {q s r±_) 2 , 

{qv/q s ) 3 Mqsr±) 

\n 2 (q s r ± ), 



r± <C q. 

q^ 1 

r 



< r± < q h 



(7) 



In this case, the In regime typically sets in only at 
extremely large length scales r± ^> q^ 1 with qh = 
q s exp[— (qy/q s ) 3 ] so that we can restrict attention to the 
In regime. At first sight, this seems to suggest the pres- 
ence of quasi-long-range translational order in this case. 
While this is formally true, the corresponding singulari- 
ties in the structure factor are extremely weak due to the 
large prefactor (qy/q s ) 3 and for practical purposes, there 
is only short-range translational order with correlation 
length ^ ans = (qs/qvY^qy 1 - As expected, the degree 
of translational order is significantly reduced compared 
to the weak-disorder regime. The degree of orientational 
order is also reduced, since g(r±) ~ (qy / k) 2 (qy / q s ) 
for r± 3> q<T . Nevertheless, £™ lcnt = 00 as long as 
q v <C k 2 q s . Otherwise, for extremely strong disorder, 

^orient ^ k 2 j q\ . 

Results for the correlation functions in the parallel di- 
rection in the asymptotic regimes can be easily obtained 
from those in the perpendicular direction by the replace- 
ment r± — > (q a /k)r«. The small factor (q s /k) <C 1 reflects 
the fact that the system is more rigid in the parallel than 
in the perpendicular direction. It is interesting to note 
that in the absence of an alignment force (B± = 0) ori- 
enting the stripes along a preferred direction, disorder de- 
stroys translational order and leads to quasi-long-range 
orientational order, even when neglecting dislocations. 

So far the density modulation has been treated assum- 
ing the absence of any topological defects. We now show 
that disorder always induces a finite density of disloca- 
tions and estimate the transverse scale L^_ lsloc on which 
dislocations become relevant and which limits the range 
of validity of the above results. The occurrence of dis- 
locations depends on a competition between the elastic 
energy cost E se \f to create a dislocation and the energy 
gain E pm in the disorder potential. The energy gain E p - m 
originates from the reduced rigidity of the charge density 
profile in the presence of dislocations. The scale i^isioc 
is identified as the size of a finite system where the free 
energy E sc u + E p - m of a single dislocation becomes nega- 
tive. For a strictly smectic system (B± — 0), the elastic 
energy cost -E S cif of a dislocation is known to be finite. 



3 



A finite tilt modulus B± leads to a logarithmic increase 
in the elastic energy with system size L± for L± ^> q^ 1 . 
The energy of the dislocation in the disorder potential 
is a random variable that depends on the location of the 
dislocation. The pinning energy E p - m is given by the max- 
imal energy that the dislocation can gain from the disor- 
der potential. Assuming a Gaussian distribution for the 
disorder energy, this energy can be computed using the 
statistics of extreme values. In the weak- disorder regime, 
we eventually find that the creation of dislocations be- 
comes favorable on transverse scales 



L 



disloc 



q v 1 exp[c\n 2/3 (q s /q v )} 



(8) 



(with c a constant of order unity) correspond- 
ing to an areal dislocation density of ridisioc ~ 
(L^ lsloc ) _2 (i3_L/i?ll) 1 / 2 . In the strong- disorder limit, we 
find 



' disloc 



- 1 min[l, k 2 q s /q v ] 



(9) 



1/2 _ 

within a finite range of disorder strengths, 

disloc 



V « Ay « (k/q s ) w / 3 A c v or q s 3 < q y < k 2 q s , L 



and a dislocation density of ridisioc ~ {L±_ s ) {K/B\\) 
Thus, 
A v 

saturates at g" 1 , before it shrinks further, coinciding 
with £°™ nt . 

The present analysis of regimes neglects the influence 
of the anharmonic elastic terms, although they are in- 
cluded at one-loop level in the flow equations (||). On 
scales q± <C qy these terms lead to an effective reduction 
of Bu and an increase of K, making the system appear 
more isotropic. The integration of the full flow equations 
p3[ modifies the correlation functions only slightly in the 
strong-disorder limit. In addition, q s is reduced due to 
the renormalization of K. 

Our results therefore lead to the following overall pic- 
ture, as summarized in Fig. 1. In the clean limit we 
have a sequence of transverse length scales g" 1 < q^ 1 ~ 



£trans ^ £ 



disloc 



(see Fig. 



On scales L± <C q v 



the system looks like a smectic with pronounced trans- 
lational order (/ <C a 2 ) and no free dislocations. On 
scales qy 1 < Li <C L^ 100 , there are still no free dislo- 
cations. However, f ^> a 2 and / increases logarithmi- 
cally, displaying quasi-long-range order, whence we call 
this regime quasi- smectic. For L± 3> L± dislocations 
appear and destroy translational order. Neglecting dis- 
locations, one finds g <C 1 (i.e. £° rlont = oo) on all scales 
in the weak-disorder regime, which indicates the possibil- 
ity to have good orientational (nematic) order even for 
Li > L 



disloc 



In the strong-disorder regime, the sequence of scales 
is ^' ans < q^ 1 < Lf sloc < q s . Good translational or- 
der survives only for Lj_ <C ^ ans - Beyond this scale, 
translational order is rapidly lost as / increases as a 
power law with / ^> 1. On the other hand, orienta- 
tional order is still retained, g <C 1. Hence, we cross 
over from a smectic regime for L± <C ^' ans to a ne- 
matic regime for L± ^> ^ ans . Dislocations are present 



on scales L± > Lf s}oc . For < Ay < (fc/9s) 10/3 A^, 
^orient — oq m the absence of dislocations and nematic 
order may survive beyond L^isioc even m their presence. 
For Ay > (fc/t7 s ) 10 / 3 A^, ^ riont ~ Lf sloc and most likely 
nematic order is destroyed for L± 3> £° rient , leading to 
an isotropic regime. 

For typical experimental samples Q , we estimate q s ~ 
qv ~ O.Olfc, indicating that they are located in the 
crossover region between the weak- and strong-disorder 
regimes (cf. Fig. 1). While we cannot firmly deduce 
the presence of orientational order on the macroscopic 
scales probed by dc transport measurements, there is 
a good chance that orientational order persists even for 
Ljl > Lf sloc , since Q ieat = oo for Ay < (fc/g s ) 10 / 3 A^ 
in the absence of dislocations. The regimes identified in 
this paper at smaller lengths may be accessible to other 
experimental probes such as surface acoustic waves. 
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